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were performed at Lawrence Livermore National 
Lab, using a Cray YMP-C90 for the Fourier ac- 
celeration and a Cray-2 to test all the other al- 
gorithms. We do not present timings of the pro- 
grams because we have not fully optimized our 
codes. 

For each algorithm we found the maximum 
value of the step size parameter a which could be 
used without destabilizing the algorithm. Figure 
1 shows div A as a function of iteration number. 
This was similiar to all the runs we did, except 
sometime the DAUB2 acceleration was a bit bet- 
ter than the unaccelerated algorithm. The prac- 
tical performance of each algorithm is a function 
of time slice. For example, on one time slice we 
tested, the DAUB6 gauge fixing algorithm only 
reduced the number of iterations by a factor of 
1.2 over the unaccelerated algorithm; for compar- 
ison, on the same time slice, Fourier acceleration 
reduced the number of iterations by a factor of 
3.4. For the wavelet transforms we tested, the ac- 
celeration improves as the wavelet gets smoother. 
The step size parameter for the unaccelerated al- 




terations 

Figure 1. Comparison of various transform accel- 
erated gauge fixing algorithms. 



gorithm (a = 0.19) was an order of magnitude 
larger than for the DAUB6 wavelet algorithm ( 
a = 0.015). This could be caused by the pre- 
conditioner we used having a destabilizing large 
step size for the biggest length scale. We also 



tried using a preconditioner which increased by 
a factor of [3 between length scales, where [3 lies 
between 1 and 2. By decreasing [3 from 2 and in- 
creasing a in an ad-hoc way, we were able to im- 
prove the performance of the wavelet algorithm. 
This suggests that there are gains to be made by 
systematically tuning the parameters of a more 
sophisticated preconditioner. 

6. CONCLUSIONS 

Although for the 8 4 lattices of this study, 
Fourier acceleration is superior to wavelet accel- 
eration, we have indicated the areas where the 
wavelet algorithms could be improved. A clearer 
picture will emerge after we have studied the 
volume dependence of all the gauge fixing algo- 
rithms. 
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Daubechies [3]. All the DAUB wavelets have com- 
pact support and are characterised by the vanish- 
ing of the moments of tp(x) 

As n increases the wavelets get bigger and 
"smoother" . The DAUB2 wavelet transform is 
the Haar transform. We remark that wavelets 
can be used to characterise smoothness spaces [3] 
(e.g. Sobolev spaces), and this could be useful 
for the development of fast algorithms. The folk 
lore in Fourier acceleration [1] is that the prac- 
tical speed of a QCD algorithm depends on how 
"smooth" the underlying gauge configuration is. 

Discrete wavelet transforms exist, with fast al- 
gorithms and published codes [6]. A Daubechies 
wavelet transform will transform a vector of 
length N in O(N) operations. We used tensor 
products of one dimensional wavelet transforms 
for the three dimensional transform. More so- 
phisticated multidimensional wavelet transforms 
exist [3], and look potentially useful. 

3. COULOMB GAUGE FIXING 

The gauge transformation G(x) which fixes a 
gauge configuration U^x) into Coulomb gauge is 
found by maximising the following function 

F = ^rw EE trace (^w + u ^) v ) ( 5 ) 

c x k = l 

where the gauge transformed configuration is 

U 3 k {x) = G(x)U k (x)G\x + k) (6) 

The algorithm used is a steepest descent method 
[1]. The gauge transformation for a single itera- 
tion is 

G(x) = exp(iw a (x)T a a) (7) 

where T a are the generators of SU{2>), a is a step 
size parameter, and the function w a {x) is defined 

by 

3 

w a {x) = -i^ trace T a A_ k (U k (x) - Ul(x)) (8) 

k = l 



The covariant derivative A_^ is defined by 

A_ k U k (x) = U k (x - k) - U k (x) (9) 

We monitor the convergence of the algorithm, 
by calculating the lattice analog of div A. 

divA = ^E^E^w 2 ( 10 ) 

X_ a 

4. WAVELET ACCELERATION 

The wavelet preconditioned gauge transform is 
G{x)jj = exp(_D~ 1 PD r iw a (x)T a a) (11) 

where D r is a wavelet transform labelled by r 
and P is a diagonal preconditioner matrix. This 
should be compared to the Fourier accelerated 
gauge transform [1] 

2 

G(x) FFT = eMF~ lE ^F iw a (x)T a a) (12) 

where F is the fast Fourier transform. 

In this exploratory study, we used the simplest 
form of preconditioner for the wavelet accelerated 
gauge transformation [7] . If the dimension of the 
space is 2 n , where n is an integer, the precondi- 
tioner matrix P in one dimension is defined as 

Pa = 

V(i) = 2 n ~ q ,2«" 1 <i<2« ,q = l---n 
V(l) = 2"- 1 (13) 

This corresponds to weighting the smallest scale 
structure with the smallest weight. The next 
biggest structure is weighted with twice the 
weight of the smallest structure, and so on, up 
the scales in the problem. 

Our wavelet transform subroutines are based 
on the routines in Numerical Recipes [6], which 
stop the recursive algorithm at the second largest 
length scale. 

5. NUMERICAL RESULTS 

We compared the performance of the different 
gauge fixing algorithms on 8 4 , (3 = 5.7, quenched 
SU{2>) gauge configurations, generated by the 
pseudo heat bath algorithm. The calculations 
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We introduce an acceleration algorithm for coulomb gauge fixing, using the compactly supported wavelets 
introduced by Daubechies. The algorithm is similar to Fourier acceleration. Our provisional numerical results 
for SU (3) on 8 4 lattices show that the acceleration based on the DAUB6 transform can reduce the number of 
iterations by a factor up to 3 over the unaccelerated algorithm. The reduction in iterations for Fourier acceleration 
is approximately a factor of 7. 



1. INTRODUCTION 

The performances of many lattice gauge the- 
ory algorithms are hampered by critical slowing 
down. The maximum stable step size of an algo- 
rithm is governed by the smallest length scale in 
the problem. Unfortunately, using the step size 
governed by the small scale structure of the solu- 
tion evolves the large scale structure very slowly 
to the final answer. In principle it is obvious how 
to improve the algorithm: use a length scale de- 
pendent step size. It is the practice which is hard. 

A technique called Fourier acceleration has 
been used to speed up the convergence of gauge 
fixing and the calculation of quark propagators in 
QCD [1]. In this method, a length scale depen- 
dent step size is implemented, using a fast Fourier 
transform to allow the step size to be a function 
of momentum. The functional form of the step 
size is chosen by studying the free field abelian 
theory. 

However, an individual configuration from a 
simulation is not translationally invariant, nor 
"smooth" , so it is not clear that the theoretical 
gains of Fourier acceleration will be achieved in 
practice. The initial work of the Cornell group, 
using Fourier accelerated algorithms to simulate 
QCD on 8 4 lattices, showed big savings in com- 
puter time. However a later study by Hockney [2] 
on the calculation of quark propagators on lattice 
volumes between 8 4 and 16 4 failed to show the ex- 
pected and crucial volume dependence of Fourier 
acceleration. In this paper, we introduce an ac- 



celeration method, similar in spirit to Fourier ac- 
celeration, based on the wavelet transform. 

2. INTRODUCTION TO WAVELETS 

Wavelets are a new field in harmonic analy- 
sis [3,4], which provide a synthesis of ideas from 
various fields including rigorous renormalisation 
group studies [5]. We shall only quote the rele- 
vant results for algorithms. 

Wavelets are a set of basis functions for square 
integrable functions, L 2 (R), which unlike the ex- 
ponential factors of the Fourier transform are si- 
multaneously localised in position and momen- 
tum space. A function f(x) £ L 2 (R) has the 
following expansion in terms of the wavelet basis 
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The ipij (x) are generated by translations and di- 
lations of a single function tp(x) 



iP ij (x) = 2H(2 i x-j) 



(2) 



The wavelet basis functions ipij{x) are orthogonal 
over different length scales, which motivates their 
use in fighting critical slowing down. 
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(3) 



The practical utility of the almost magical equa- 
tion 3 rests on the choice of the function ip{x). 
In this work we used the wavelet family due to 



